A Creation Operator for Spinons in One Dimension 
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We propose a definition for a creation operator for the spinon, the fractional statistics elementary 
excitation of the Haldane-Shastry model, and give numerical and analytical evidence that our oper- 
ator creates a single spinon with nearly unit amplitude in the ISE model. We then discuss how the 
operator is useful in more general contexts such as studying the underlying spinons of other spin 
chain models, like the XXX- and XY-model, and of the one dimensional Hubbard model. 



One of the deepest open problems in condensed matter 
physics today is that of understanding fractional statis- 
tics, that is understanding field theories whose underlying 
excitations do not obey either Bose or Fermi statistics. 
Such possibilities (fractional statistics) for point parti- 
cles are forbidden in three dimensions; however, in, two 
dimensions, braid statistics can occur and anyonsEl are 
allowed, ^odiile, in one dimension, there exist parafermion 
operators! 2 ] which also do not obey Fermi or Bose statis- 
tics. Such fractional statistics objects may be important 
for describing the low energy excitations of condensed 
matter systems if these systems are categorized by fixed 
points whose effective, low energy Hamiltonians are one 
or two dimensional. 

The-jHeisenberg model with inverse squared exchange 
(ISE)HH (H = ^(f suitt^)"^ ■ Sj) is believed to 
provide a model whose excitations are jjaqmpletely non- 
interacting fractional statistics particlescm referred to as 
spinons. The statistics of these excitations are most nat- 
urally discussed in-jterms of an alternative definition, pro- 
posed by HaldaneQ, in which a particle's statistics is de- 
fined in terms of its effect on the size of Hilbert space for 
other particles. For example, on the lattice, each fermion 
added to a system decreases the number of states in the 
Hilbert space for other fermions of the same type by one, 
while each added boson leaves the size of the available 
Hilbert space unchanged. Fractional exclusion statistics 
occurs for operators which reduce the Hilbert space by 
a non-integer amount, i.e. at a fraction of the rate for 
fermions. This definition has been shown to agree with 
the definition of statistics through particle interchange in 
the case of the quasiparticles of the Laughlin state, (when 
the cut-offs in the theory are treated carefully )LTeI. In 
that case, the change in the phase of the wavefunction 
resulting from quasihole interchange is the same fraction 
of 7T as the rate of Hilbert space reduction due to the 
creation of quasiholes, relative to the rate of reduction 
due to the creation of fermions. 

For the exclusion approach, lattice models are more 
natural since the number of states in the Hilbert space is 
easily defined, and, indeed, Haldane introduced his def- 
inition partially motivated by results for the ISE spin 
chain model: the spectrum of the ISE spin chain model 
can be described with a set of non-interacting single par- 



ticle energies, and a set of rules for the occupation of the 
single particle levels which embody the Yangian symme- 
try of the modeH. The effect of these rules is that the 
creation of spinons eliminates allowed spinon states (of 
any polarization) from the Hilbert space half as fast as 
the creation of fermions eliminates allowed fermion single 
particle states from a Hilbert space of fermionic statesEj. 
The spinons are thus half-fermions or semions and are 
completely free. This property makes the ISE model 
somewhat more natural for studying exclusion statistics 
than the Laughlin states since the quaslholes are n °t ideal 
fractional exclusion statistics particlesO. 

In light of this, it would be very useful to have some 
method of treating the underlying spinons in that and 
related models beyond the first quantized formalism; one 
would like to introduce spinon creation and annihilation 
operators. The main purpose of this work is to introduce 
an operator which has essentially unit overlap with the 
spinon creation operator in the ISE model, and which, 
we conjecture, can be interpreted as a spinon creation 
operator much more broadly than this. 

The operator we now propose was motivated chiefly 
by the beharaor of the large U wavefunction of the Hub- 
bard modeled, which factorizes into spin and charge parts 
which are linked only by a backflow condition!! 2 ]. Except 
for the backflow condition, the charge part is described by 
free spinless fermion wavefunction and the spin part by 
a wavefunction of the one dimensional Heisenberg model 
(H = J2i Si • Si+i) where the role of the sites is played 
by the electrons. The creation operator for an up spin 
electron in that model involves the creation of a spinless 
fermion in the charge wavefunction, the backflow and 
the insertion of a site with an up spin into the Heisen- 
berg model wavefunction. Likewise, the bosonization ap- 
proach to the Hubbard modelll exhibits a similar spin 
charge decoupling, with the electron creation operator 
effectively factorizing into a product of exponentials in 
spin and charge density eigenexcitations. This suggests 
the identification of the spin insertion with the spin boson 
part of the electron operator. Further, the spin operator 
is identical to the primary field of the c = 1 conformal 
field theory studied irO, and the Fourier models of this 
operator are also recognizable as a semionic parafermion. 
Thus, it has the essential properties expected of a spinon 
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operator, being both semionic-.and providing a natural 
representation of the YangianliJ. 

We therefore conjecture that the local up-spinon cre- 
ation operator might be given approximately by an op- 
erator, 0^(x), which inserts an additional up-spin into 
a one-dimensional spin chain after site x but before site 
x + 1. Peculiarly this operator, in the process of cre- 
ating a spinon, would have to change the length of the 
spin chain it acted on from N to N+1. Some such 
an exotic effect is clearly necessary in defining a single 
spinon creation operator because of the fractional ex- 
clusion statistics of the spinon; an operator which cre- 
ates a single spinon must involve the rearrangement of 
the Hilbert space different from what occurs for fermions 
or bosons and hence the need to change the number of 
sites in the model in order to create a single spinon. We 
have discussed this feature of exclusion statistics,-before 
in a slightly different context involving spin chainsLLa, and 
it also occurs in the case of Laughlin wavefunctions re- 
stricted, for example, to the sphere where the creation 
operator involves an addition of allowed states to the 
(appropriately cut-off) Hilbert space. 

If we accept this unusual property as an unavoidable 
complication for fractional statistics operators, the evi- 
dence is there to support the proposal that the spin inser- 
tion creates a single spinon with the same polarization. 
Let us first study this conjecture in the context of the 
ISE model. In complete analogy with the insertion of two 
spins in a singlett3, we can construct the state obtained 
by inserting a single spin. To fix notation: the unnor- 
malized groundstate wavefunction on N (even) sites of 
the ISE model in a basis labeled by the positions of the 
down-spins {m, ■ ■ ■ , tim} , M = y is given by: 
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(1) 



In a basia of local spins 
duces tol_3: 



{|<7i • • -crjv)}, crj = ±5 this re- 



isgn(eri-o-j) 



(2) 



For the even N-site ISE groundstate {z^ = Cm = 
| e -|pj n — i.,, TV}. The action of 0^(x — N) on this 
state is, in the language of eq. (Q), to add an extra spin 
o~n+i =T, and leave the other spins alone (with trivial 
generalization to other values of the insertion-site x). 
In order to answer the question, to what extent this 

state Cfl (x)]^^) simulates a true-JpSpinon state, let us 
recall what the latter looks likctaO Since the ISE- 
Hamiltonian commutes with momentum, the 1-spinon 
eigenstates (which exist only for an odd number of sites) 
have definite momentum. However these states can be co- 
herently superimposed (Fourier transformed) into states 
in which the single spinon is completely localized on a 



site, say x, which is no longer an eigenstate-pof Hise- 
Such an unnormalized state is known to bet£rEII of the 
form 

Vy +1 (ni...ni*\x) = 

N_ N_ 

f[(-rM^)f[sin%^,), (3) 



in the language of down-spin co-ordinates {n^}. In terms 
of spin co-ordinates {eri . . . cttv+i} it looks exactly like eq. 

d)E 

N+1 

*£+i= II (Zi-ztf-'ietwto-rt, (4) 
but a x is excluded from the product (since o~ x is fixed 

27T1 

to be up), and the set {zi} = Cn+i = {e N + 1 ; n 
I ... A' ■ 1}/{z~>^}. Notice that in the {crj-basis the 
quasihole-factor in (|j|) goes away. Our conjecture is 

now that \^ N ] +1 (x)) ~ 0^{x)\^ { ° ] ). Numerically it 
turns out that the two states have an excellent overlap: 

(*$ +1 (x)\oHx)\*f?) 



0.984+O(7V-2). This is unsurpris- 



ing, since 

(x)) is obtained from Ot(a;)|*^ ) > b y adi " 
abatically deforming the set of N z r 's: e 2 ^ L ' r — > e JV + 1 ' r , 
and leaving the single spin(on) at site x alone, see Fig. ||. 
Although this is an uncontrolled approximation, it seems 
justified by the high numerical overlap. 

Naively one might expect from eq. (^) there to be 
[N + 1) 1-spinon states; there is, however, a reduction, 
because only N £ 2 of these are linearly independent. This 
can be seen by reconstructing the momentum eigenstates, 
through Fourier transforming eq. (^|) in the spinon co- 
ordinate x. This transform only has support in the range 
[-|, |] (N = 4m) or [f , ^] (N = 4m + 2). One notices 
that the states thus obtained are indeed eigenstates, since 
(1) the set {^^ +1 (x)\x = 0...N} spans the 1-spinon 
subspacal3, and (2) there is only one up-spinon eigen- 
state for a fixed value of momentum. 

As stated before, spinons in eigenstates are not local- 
ized but have a momentum so that we must investigate 
to what extent our spin insertion can be connected with 
spinon creation in momentum space. Since the number 
of sites in the chains before and after the action of the 
spinon operator is different, the definition of the spinon 
creation operator of fixed momentum is not trivial and 
the Fourier transform of the creation operator requires 
some care. If we define the momentum space version of 
the creation operator by: 



0*(k 



N+1 



Qn) 



y/N 



1 N 

== y^ e i{kH+x-m)x Ux) (5) 



where k is an allowed momentum in the N+1 site model 
and q is an allowed momentum in the N site state, then, 
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the operator takes a momentum eigenstate of the iV site 
model to an momentum eigenstate of the N + 1 site 
model. (In the following we will assume periodic bound- 
ary conditions, but other boundary conditions can be 
implemented straightforwardly). In particular, if N is a 
multiple of 4 then the groundstate has zero momentum 
and q can be taken to be zero if we wish to create a single 
excited spinon. In this case, this construction allows us 
to study our approximation to the single spinon spectral 
function defined by: 

A s (k,u) = Y / \WOHk N +i)\GS)\ 2 6(E m -E -uj) (6) 



where the sum over m is over a complete set of eigenstates 
of the N + 1 site model and \GS) is the groundstate of 
the iV site model. If our conjecture were a perfect spinon 
creation operator, then for the ISE model the spectral 
function would be a (^-function of unit weight located on 
the 1 spinon dispersion relation e s (fc) = ^((tt/2) 2 — k 2 ). 

The actual result is remarkably close to this as shown 
in Figure ||. For the finite systems that we studied (up to 
24 sites), for the ISE model, at least 99.4% of the weight 
of the spin insertion seems to lie on the single spinon. 
However the weight in the spectral function for O* is not 
unity (as it would be for ct in the case of free bosons 
or fermions), nor is it close to it, since we don't know a 
priori how to normalize our spinon operator such that it 
creates a normalized state. However, for the ISE, using 
the aforementioned almost-equivalence of spin-inserted 
wavefunction and the localized spinon wavefunction, we 
can get an idea of what this normalization should be like. 
That is, we will find the normalization appropriate if the 
spin-insertion put all of its weight in the onc-spinon state 
(rather than just 99%). Then we need only compute 
the contribution to the normalization from the matrix 
element in (^) , for the 1 spinon state (of momentum Q = 
k). The weight A(k,s s (k)) is given by (in an obvious 
notation): 



N 



VNC^CWA(k) = (¥$ +1 (Q = • £ l*#n(j)>e- ttj 

3=0 

N 

= E^+iWi^+iC?)}^*-^, (7) 

i,j=a 

where C<°) = $ N m 2 -m(n-i) 

and CW = M \ N m 2 -m(n~i) are the norma ii zat ion of 
respectively states (eq. |l|) and ^^ +1 (x) (eq. ||). 
The relevant matrix element in the RHS of this equation, 
involving sums over locations of down-spins, can be com- 
puted, after one realizes that these sums can be replaced 
by integrals of the same expression, due to the polyno- 
mial nature of the summandHjIZl. Methods to compute 
integrals of this kind have been developed in the con- 
text of the Calogero-Sutherland modeH and random ma- 
trices, using for instance Jack-polynomials, and matrix- 



model correlators EES. We will follow instead Suther- 
land's recipea,E3. There we find, after a little bit of alge- 



bra that (M = f 



M 



i=i 



^) = ±(c* S (^7r)-COs(0)) 

u(9) {p cos 9p cos Oq + q sinp6> sin q9} 

^2irp 

= 1 | cos (nigi) S pq - ^S M+X - ~<W-i } (*+f ) 



and p, q = - . . . M - ^ 



(8) 



If we set x = cos (^tt) the ensuing (tridiagonal!) de- 
terminant observes the recursion relation: 



d n +i(x) = {2n + l)d n (x) - n 2 d n ^i(x) 



(9) 



with initial conditions do = 1, d\ — x. This is exactly 
the definition of the Legendre Polynomials: oIm{x) = 
M\Pm{x) (with normalization Pjvf(l) y=\ 1). Then the 
Fourier transform in eq. (0) is given by 



A(k) 



(N + l) 



2r \ ( 2(M - r) \ 9 _ 2 M 



M-r 



\/(§) 2 - fc2 ' 



(10) 



and k = ^ — (r+j). In the derivation we have chosen 
M even to prevent cumbersome notation. See the remark 
above about where the weight is located, depending on 
the parity of M. We notice that the weight diverges at 
what Haldane calls the spinon "pseudo-fermi surface" . 
The square root divergence is exactly that expected for 
spectral function of the spin part of the electron operator 
obtained in Abelian bosonization; the exponent shouldJae 
-^p + j^- = h, for K a = 1, i.e. isotropic spin exchangees. 

Although almost all of the effect of the operator defined 
by eq. is to create a single spinon in an ISE models, 
there is some amplitude in the spectral function off the 
single spinon energy, so that our operator can not be 
taken to be identically the spinon creation operator. The 
connection between the two is rather like that between 
the bare electron operator and the quasiparticle operator 
in Fermi liquid theory: in addition to the expected 5- 
function there are very small "incoherent" contributions 
coming from the fact that our creation operator has finite 
weight to create more than one spinon. 

In any case, the operator will be extremely useful if we 
can demonstrate that it has large finite overlap with the 
spinon creation operator in a broader context than acting 
on the ground state of the ISE model. For example, the 
Heisenberg model is in the same universality class as the 
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ISE model and can be thought of as a model of nearly 
free spinons (they have a marginally irrelevant interac- 
tion), and so we have displayed the spectral function of 
the spinon operator acting on the groundstate of that 
model together with the ISE spectral function in Fig- 
ure ||. For both models there is very little probability to 
excite states with momentum, inside the "pseudo-fermi 
surface" , weight is down there by as much as 3 orders of 
magnitude, and, we see that, for the Heisenberg model, 
the weight is predominantly distributed on a 1-parameter 
family of states, which carries at least 98 % of the weight. 
Notice that the single spinon feature in the spectral func- 
tion exists only over half the Brillouin zone as expected 
from the Bethe Ansatz solution of .the Heisenberg model 
and the solution of the ISE modelD. This is also consis- 
tent with the fact that the two species of spinon (up and 
down spin) should obey mutual | statistical so that the 
creation a spinon of either type reduces the Hilbert space 
of the spinon of the other type by ^ . This is realized here 
since adding one site and one spinon to the model adds 
only i of a state for other spinons, because only one-half 
of the Brillouin zone is accessible. The interchange of 
two spinon insertion operators acting at different points, 
conversely does not yield a phase of ±7r/2; rather the 
leftmost spinon is shifted one site. However, in the low 
energy limit, the spinon lives near momentum ±7r/2 and 
the position shift mimics semionic statistics. 

Both the ISE and Heisenberg models are exactly sol- 
uble and therefore well understood, however, our defini- 
tion can be extended to essentially any one dimensional 
spin chain model, including anisotropic spin chain model 
such as the XY model (H = £\ SfSf +1 + SfSf +1 ). This 
model can be mapped onto a hard-core boson problem 
and then onto free fermions via a Jordan- Wigner trans- 
formation, so it is readily soluble, however, the connec- 
tion between that solution and other viewpoints is still 
unclear. In particular, it should be possible to under- 
stand the XY model as model of interacting spinons. If 
our spinon creation operator can really be taken as such 
then we can study the underlying, interacting spinons 
of XY model and elucidate this connection — a possibil- 
ity that would not exist were we restricted to the first 
quantized understanding of the spinons available for the 
ISE model. For example, we have computed the spinon 
spectral function as before for this model with the results 
shown in Figure |^. 

They can be compared to the theoretical explanation 
based on the identification of the chiral parafermion, 
exp(i^9^ or L (x)), with the spinon. In that case, the 

Luttinger liquid hypothesis!! for the XXZ chains pre- 
dicts that the asymptotic correlations of the spinon in 
the XY model should take the form 



-^spinon (*E : ^) 



^iirx /2 



computed from thisc-3 and it follows that: (1) the spinon 
spectral function has no deZta-functions but rather power 
law singularities at lu — ±vk, with k measured from 
±tt/2, of the form: \u - vk\~ 15 / 16 and \u> + vk\' 7 / 16 and 
(2) the integral of the spectral function over frequency: 
A(k) (which would simply be the 1 — for fermions or 
bosons) diverges like k~» as ±7r/2 is approached from 
either side. The ratio of the prefactors of the divergences 
5.03. This agrees with our numerical re- 



sin(77r/16) 
sin(7r/16) 



(X - Vt) 9 / 16 (x + vt)T6 

left moving piece + non — universal (11) 



The universal properties of the spectral function can be 
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suits for this ratio for systems of large size (about 2000 
sites for A{k) and 100 for A spinon (fc, u)). The power law 
rather than <5-function divergence in Aspi non is quite ob- 
vious when we compare Fig. | and §: ie. O^(Q) has 
significant weight on a large number of states beyond the 
lowest-energy single particle-hole state. One apparent ex- 
ception to this statement is spin-insertion with Q ~ 
in that case the larger part of the weight remains on a 
single particle-hole state, even as N — ► oo. However, 
the relative weight on this class of free Jordan- Wigner 
fermion states decays as [|f — Ql^\ r ]~ r ' with r\ ~ 0.1 and 
vanishes (slowly) when we take N — > oo before \Q\ — > ?. 
Thus our spinon construction is in accord with all of the 
expected properties for a spinon creation operator in the 
XY model. 

An analog of our spinon operator can also be con- 
structed for "spin" chains not based on SU(2), such as 
the SU(3) generalization of the nearest neighbor Heisen- 
berg model in which the "spin" at .each site is in the fun- 
damental representation of SU(3)c]i (This is the model 
that would naturally be obtained from an SU (3) gener- 
alization of the large U Hubbard model at ^ filling). The 
results for the spectral function are shown in Fig. 

For the ISE and Heisenberg models, we knew that we 
were dealing with free or weakly interacting spinons, how- 
ever, the existence of such spinons in the SU (3) model, 
while plausible has not been shown. The spinon spectral 
function of Fig. [| demonstrates that that model is also 
one of nearly free spinons. Our spinon creation operator 
construction is thus very useful in this context, yielding 
a powerful, qualitative result quite simply. Moreover, 
the spectral function we have obtained reveals that the 
spinons of the SU(3) model are not semions but rather 
obey mutual ^ statistics; that is the creation of one of any 
of the three spinon species reduces the Hilbert space for 
any of the three spinon species by | of a state. This fol- 
lows from the fact that the "spinon" spectral function has 
support only over | (in general for SU(p) spin mod- 
els) of the Brillouin zone. These results aijedn agreement 
with findings for the ISE SU(N) modelscl and demon- 
strate the the nearest neighbor model is almost certainly 
in the same universality class. 

Our results for the spinon spectral function in the 
Heisenberg model are directly related to the work on the 
electron sp ec t ra l function of the one dimensional Hub- 
bard mode£3c3. As a result of the completeness of the 
momentum space spinon operators for states reached by 
the insertion of a single spin, and the fact that our spinon 
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creation operator creates a momentum eigenstate, the 
function C a (k, u>) ofcJis, in fact, identical to our A s {k,w). 
Their finding that the weight in C a (k,uj) is concentrated 
in half the Brillouin zone at the lowest allowed energies is 
a consequence of the fact that C a (k,Lu) contains the sin- 
gle spinon spectral function plus small corrections com- 
ing from three and higher spinon terms. As a result, 
it is dominated by the single spinon part of the spectral 
function. This occupies only half the Brillouin as a conse- 
quence of the Yangian symmetry and fractional statistics 
of the spinonsB. For the XY model, (also studied inn3) 
the spinons are interacting and even the single spinon 
contribution to the spectral function has weight over the 
entire Brillouin zone and a power law divergence as the 
threshold energy is approached from above. This be- 
hajiior for C a (k, u>), was, in fact, independently obtained 
irM 

In summary, we have proposed a spinon creation op- 
erator for one dimensional spin models and their gener- 
alizations. We have shown that the states created by the 
proposed operator have excellent, although not perfect, 
overlap with the actual one spinon eigenstates of the ISE 
model. In addition, for finite size Heisenberg models, 
nearly all of the operator's weight when acting on the 
groundstate goes into creating an the eigenstate with the 
lowest possible energy for a given momentum — consistent 
with what expects for the single spinon creation operator 
given that the Heisenberg model is a model of spinons 
with marginally irrelevant interactions. For the XY 
model, the operator acting on the groundstate creates 
states with a broad distribution of energies, but the de- 
tailed properties of the distribution are those expected 
for a spinon creation operator (if the XY model is re- 
garded as a Luttinger liquid model of strongly interact- 
ing semionic spinons). For an SU(3) generalization of 
the Heisenberg model, we find results consistent with a 
straightforward generalization of the SU(2) spinon cre- 
ation operator acting as a spinon creation operator for a 
system of weakly interacting "spinons" obeying Haldane 
type ^ mutual exclusion statistics. This result is con- 
sistent with expectations but, to our knowledge, has not 
been demonstrated by any other means. 

Together, these results indicate that the operator pro- 
posed is a valid spinon creation operator, and should be 
quite useful in the study of one dimensional spin models 
and their generalizations. 

One of us (S.P.S) acknowledges support from Depart- 
ment of Energy grant DOE DE-FG02-90ER40542. 
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FIG. 1. Graphical illustration of the deformation of the 
spin insertion to a spinon singlet insertion. The gray area 
is the groundstate wavefunction, or an adiabatic deformation 
thereof. 




FIG. 2. The spectral weight of the spin-insertion operator for (a) ISE and (b) Heisenberg ID spin-chains, on 16 sites. The 
discrete data has been smoothed by convolving with a Lorentzian. Notice that the z-axis has a logarithmic-scale. Otherwise, 
on a linear scale, the multi-spinon contributions would not be visible! The solid line on the back wall, A(k), is the ^-integrated 
spectral weight, or equal time spin- insert ion correlator. The dashed line is the contribution from single spinon states (the height 
of the large peaks) . On the top we show a density-plot from which it is clear that virtually all the weight lies on the 1-spinon 
dispersion curve — e B (k) = ) 2 — k 2 ) for ISE and 7rcos(fc) for the Heisenberg modelu. 
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FIG. 3. The spectral weight of the spin-insertion operator 
for the XY-model on 100 sites, allowing up to 3 particle hole 
excitations. The data has been smoothed by convolving with 
with a Lee-filter, and plotted on a log-linear scale. Notice that 
some finite weight seems to be present on the groundstates 
with momenta — 5 < Q < 5; these contributions, however, 
vanish in the thermodynamic limit, except for the lowest en- 
ergy states close to ±-|. The latter have measure 0, thus 
restoring results from bosonization. The curve on top is the 
1 particle-hole dispersion, other symbols are as in fig. 0. 




FIG. 4. The spectral weight of the spin-insertion operator 
for the SU(3) Heisenberg model on 9 sites. The data has been 
smoothed by convolving with a Lorentzian and plotted on a 
log-linear scale. The symbols are as in fig. ^. Notice that 
in addition to the 1-spinon (^'-representation) there is also 
some weight on states with 2 anti-spinons ('3'rfepresentation). 
This feature remains presumably as N — > oc£3. 
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